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ABSTRACT 

We have analysed 53 square degrees of i^c-band imaging data from the Red-Sequence Cluster Survey 
(RCS), and measured the excess correlations in the shapes of galaxies on scales out to ~ 1.5 degrees. 
We separate the signal into an "E"- (lensing) and "B"-mode (systematics), which allows us to study the 
contribution from residual systematics and intrinsic alignments. On scales larger than 10 arcminutes, 
we find no "B"-mode, suggesting that the signal at those scales is solely caused by gravitational lensing. 
On smaller scales we find a small, but significant "B"-mode. This signal is also present when we select 
a sample of bright (20 < Rc < 22) galaxies. These galaxies are rather insensitive to observational 
distortions, and we therefore conclude that the observed "B"-mode is likely to be caused by intrinsic 
alignments. To minimize the effect of intrinsic alignments, we limit the cosmic shear analysis to galaxies 
with 22 < Rc < 24. 

We derive joint constraints on flm and erg, by marginalizing over F, Q\ and the source redshift 
distribution, using different priors. Marginalizing over F and JIa, and using a flat prior for the source 
redshift distribution, yields a conservative constraint of cts = 0.4:5t.Q'^2^^'^^ (95% confidence). A 
better constraint is derived when we use Gaussian priors for F (from the 2dF survey) and il™ + ^Ia (from 
CMB), and the source redshift distribution. For this choice of priors, we find as = (0.46t°:^^)r!-''-52 
(95% confidence). We also investigated whether the RCS data can be used to constrain F. Using our set 
of Gaussian priors, we find that we can only place a lower bound on F for which we find F > O.I + 0.I6r2m 
(95% confidence). 

Comparison of the RCS results with three other recent cosmic shear measurements shows excellent 
agreement. The current weak lensing results are also in good agreement with CMB measurements, 
when we allow the reionization optical depth r and the spectral index Us to vary. We present a simple 
demonstration of how the weak lensing results can be used as a prior in the parameter estimation from 
CMB measurements to derive constraints on the reionization optical depth r. 

Subject headings: cosmology: observations — dark matter — gravitational lensing 



1. INTRODUCTION 

The measurement of the coherent distortions of the im- 
ages of faint galaxies caused by weak lensing by intervening 
large scale structures provides a direct way to study the 
statistical properties of the growth of structures in the uni- 
verse. Compared to many other methods, such as redshift 
surveys (e.g., Strauss & Willick 1995; Peacock et al. 2001; 
Percival et al. 2001), weak lensing probes the dark matter 
distribution directly, regardless of the light distribution. 
In addition, it provides measurements on scales from the 
quasi-linear to the non-linear regime, where comparisons 
between observations and predictions are still limited. 

The distortion in the images of distant galaxies is small, 
and consequently it is not measured easily. Only recently 
it has become possible to measure the signal and the 
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number of published measurements is increasing rapidly 
(e.g.. Bacon et al. 2000, 2002; Hoekstra et al. 2002a; 
Kaiser, Wilson, & Luppino 2000; Maoh et al. 2000; Re- 
fregier, Rhodes, & Groth 2002; van Waerbeke 2000, 2001a, 
2002; Wittman et al. 2000). 

Weak lensing observations can be used for a wide range 
of studies. For instance, Hoekstra, Yee, & Gladders (2001) 
presented a first measurement of a combination of the 
bias parameter and the galaxy-mass cross-correlation co- 
efficient as a function of scale. The first direct measure- 
ment of the bias parameter and the galaxy-mass cross- 
correlation coefficient was presented by Hoekstra et al. 
(2002b). These studies provide unique constraints on mod- 
els of galaxy formation. Another important application is 
to constrain cosmological parameters from the measure- 
ment of the two-point statistics of the dark matter distri- 
bution. 

The amplitude of the lensing signal depends on many 
parameters, of which some are degenerate. Some of the 
degeneracies can be overcome using measurements that 
probe the dark matter power spectrum at different red- 
shifts, or by comparison of the results in the linear regime 
to the non-linear scales. Such subtle differences require 
very accurate measurements, which are currently not avail- 
able. Planned large surveys, such as the CFHT Legacy 
Survey, however, will be a major advance in this area. 
Also higher order statistics, such as the skewness, are po- 
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tentially powerful tools to break some of the degeneracies 
(e.g., Schneider ct al. 1998). 

Better constraints can be obtained by combining the 
weak lensing measurements with results from other tech- 
niques. In particular the combination with cosmic mi- 
crowave background (CMB) measurements is useful, be- 
cause the physics of both lensing and CMB are well un- 
derstood. Also some results from galaxy redshift surveys 
can be used, provided they are made in the linear regime, 
and under the assumption that the bias is constant with 
scale. 

In this paper we present the results from our analysis of 
~ 53 deg^ of iJc-band data from the Red-Sequence Clus- 
ter Sequence (RCS) (e.g., Yee & Gladders 2001), which 
is a 90 deg^ galaxy cluster survey designed to provide a 
large sample of optically selected clusters of galaxies with 
redshifts 0.1 < z < 1.4. 

Hoekstra et al. (2002a) used a subset of the current 
data to study the cosmic shear. Compared to other cos- 
mic shear studies, the RCS data arc shallow, and conse- 
quently the signal at a given scale is much lower, as is the 
signal-to-noise ratio. However, measuring the weak lens- 
ing signal from a shallow survey has several advantages. 
Down to a limiting magnitude of Rc ~ 24 star-galaxy sep- 
aration works well (see Fig. 1 from Hoekstra et al 2002a). 
In deeper surveys many sources have sizes comparable to 
the size of the PSF, and applying size cuts may change the 
redshift distribution of the sources in a systematic way. In 
addition, down to Rc ^ 24 the redshift distribution of the 
sources is fairly well determined. In order to relate the 
observed cosmic shear signal to cosmological parameters, 
a good understanding of the source redshift distribution is 
crucial. 

One worry is the effect of intrinsic alignments of the 
source galaxies, which introduces an additional signal (e.g.. 
Heavens et al. 2000; Catelan et al. 2001; Crittenden et 
al. 2001; Mackey et al. 2001). The amplitude of the ef- 
fect is not well determined, but it is clear that it is more 
important for shallower surveys. However, the predictions 
indicate that for a median redshift of z = 0.6 (which is sim- 
ilar to our sample of source galaxies) the signal caused by 
intrinsic alignments is still small compared to the lensing 
signal. 

The structure of the paper is as follows. In §2 we briefly 
discuss the theory of lensing by large scale structure. The 
data and measurements are presented in §3. We compare 
the measurements to model predictions in a maximum like- 
lihood analysis in §4. We summarize our results in §5. 

2. METHOD 

In this section we provide a basic description of the the- 
ory of weak lensing by large scale structure. We discuss the 
dependence on the assumed cosmology and the redshifts 
of the sources. Detailed discussions on this subject can be 
found elsewhere (e.g., Schneider et al. 1998; Bartelmann 
& Schneider 2001). 

The observable two-point statistics can be related to the 
convergence power spectrum, which is defined as 
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where w is the radial (comoving) coordinate, Wh corre- 
sponds to the horizon, a(w;) the cosmic scale factor, and 
fxiw) the comoving angular diameter distance. As shown 
by Jain & Seljak (1997) and Schneider et al. (1998) it is 
necessary to use the non-linear power spectrum in equa- 
tion (1). This power spectrum is derived from the linear 
power spectrum following the prescriptions from Peacock 
& Dodds (1996). 

W{'w) is the source-averaged ratio of angular diame- 
ter distances Dis/Dg for a redshift distribution of sources 
Pb{w): 
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Hence, it is important to know the redshift distribution 
of the sources, in order to relate the observed lensing signal 
to F^l). 

In Hoekstra et al. (2002a) we used the top-hat smoothed 
variance as a measure of the cosmic shear signal. The 
top-hat variance (7^) is related to the convergence power 
spectrum through 
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where 6 is the radius of the aperture used to compute the 
variance, and Ji is the first Bessel function of the first kind. 
Hoekstra et al. (2002a) measured the top-hat variance by 
actiially computing the excess variance in apertures. 

Recent studies, however, show that an optimal use of the 
data is to measure the shear correlation functions from the 
data directly. These correlation functions can be related 
to the various two-point statistics. In addition, this ap- 
proach allows one to split the signal into two components: 
an "E"-mode, which is curl- free, and a "B"-mode, which 
is sensitive to the curl of the shear field. Gravitational 
lensing arises from a gravitational potential, and hence it 
is expected to produce a curl-free shear field. Hence, the 
"B"-mode can be uscid to quantify the level of systematics 
involved in the measurement. 

Several sources of "B"-mode have been identified. 
Schneider et al. (2002) showed how redshift clustering 
of source galaxies can introduce a "B"-mode. However, 
in most situations the contribution to the observed sig- 
nal is negligible. Simple models describing the intrinsic 
alignments of galaxies predict a small "B"-mode (Critten- 
den et al. 2002), although the amplitude is still uncertain. 
Hence, any measured "B"-mode is dominated by residual 
systematics in the data (e.g., imperfect correction of the 
PSF anisotropy) or intrinsic alignments. The effect of in- 
trinsic alignments can be minimized by selecting galaxies 
with a broad redshift distribution. In future surveys, with 
photometric redshift information for the galaxies, the con- 
tribution from intrinsic alignments can be removed com- 
pletely by correlating the shapes of galaxies with different 
redshifts. 

The decompositions of the shear correlation function 
and the top-hat variance in E and B modes are defined 
up to a constant (Crittenden et al. 2002; Pen et al. 2002). 
The decomposition is naturally carried out by using the 
aperture mass statistic Map. Kaiser et al. (1994) intro- 
duced this statistic to measure the masses of clusters of 
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galaxies. Its usefulness for the study of cosmic shear was 
pointed out by Schneider et al. (1998). The aperture mass 
is defined as 

Map(^) = J d^cf>U{<P)K{<f>). (4) 

Provided U{(j)) is a compensated filter, i.e., J d<p(l)U{<p) = 
0, with U{(p) = for (/) > 6, the aperture mass can be 
expressed in term of the observable tangential shear 7t us- 
ing a different filter function Q{(p) (which is a function of 



(5) 



Map(0) = / d^0g(</.)7t(0). 



We use the filter function suggested by Schneider et al. 

(1998) 



(6) 



for 9 < d, and elsewhere. The corresponding Q{9) is 
given by 
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for < {}, and elsewhere. With this choice of filter func- 
tions, the variance of the aperture mass (Mfp) is related 
to the power spectrum through 
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where J4 is the fourth-order Bessel function of the first 
kind. 

A straightforward implementation of the aperture mass 
statistic is to compute Map directly from the data using 
the estimator (e.g., Schneider et al. 1998) 
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where Ns is the number of source galaxies, 7t,i is the tan- 
gential shear of the ith galaxy with respect to the center 
of the aperture. The weights Wi correspond to the in- 
verse square of th(^ uncertainty in the shape measurement 
(Hoekstra et al. 2000). 

One can use this procedure to tile the observed fields 
with apertures, and compute the excess variance. How- 
ever, this approach assumes that the data are contiguous, 
which is not the case for real data. The masking of the 
RCS data is not severe, but does complicate a measure- 
ment of (M^p) through this method. 

Instead, we will determine (M^p) from the observed el- 
lipticity correlation functions. This method has the advan- 
tage that it uses all information contained in the data, and 
that it does not depend on the geometry of the survey. Pen 
et al. (2002) were the first to apply this approach, using 
data from the VIRMOS-DESCART survey. 

The two ellipticity correlation function that are mea- 
sured are 
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where = — Xj|. 7t and 7]. are the tangential and 45 
degree rotated shear in the frame defined by the line con- 
necting the pair of galaxies. For the following, it is more 
useful to consider 



e+(^?) =eu(e)+4r(e), andC-(0)=ett(e)-4r(e). (12) 
i.e., the sum and the difference of the two observed cor- 
relation functions. As shown by Crittenden et al. (2002), 
one can derive "E" and "B"-mode correlation functions by 
integrating S,+ {0) and (,-{0) with an appropriate window 
function (see Pen et al. 2002 for an application to the 
VIRMOS-DESCART data). 

The "E" and "B"-mode aperture masses are computed 
from the ellipticity correlation functions using 
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where W^d), and W(??) are given in Crittenden et al. 
(2002). Useful analytic expressions were derived by 
Schneider et al. (2001). Both >V(i9), and W{d) vanish 
for -d > 26, so that (Affp) can be obtained directly from 
the observable ellipticity correlation functions over a finite 
interval. 

3. MEASUREMENTS 

We use the Rc-hmid imaging data from the Red- 
Sequence Cluster Survey (RCS, Gladders & Yee 2000). 
The final survey covers 90 deg^ in both Rc and z', spread 
over 22 widely separated patches of ^ 2.1 x 2.3 degrees. 
The northern half of the survey was observed using the 
CFH12k camera on the CFHT, and the data from the 
southern half were obtained using the Mosaic II camera 
on the CTIO 4m telescope. 

The integration times arc 900s for the CFHT data and 
1200s for the CTIO data. A subset of this data set was 
studied in Hoekstra et al. (2002a): 16.4 deg^ CFHT data 
and 7.6 deg^ of CTIO data. We augment the data used 
by Hoekstra et al. (2002a) by the remaining CFHT12k 
data. Much of the remaining CTIO data has relatively 
poor seeing, and hence is less useful for our weak lensing 
analysis and is not used. The seeing for the data used here 
ranges from 0'.'57 to 1'.'39 (only 11% of the data had seeing 
> 1"), with a median value of 0'.'8. Hence, the final data 
set studied here consists of a total of 53 deg^ of imaging 
data, spread over 13 patches. 

A detailed description of the data reduction and object 
analysis is described in Hoekstra et al. (2002a), to which 
we refer for technical details. Here we present a short de- 
scription of the various steps in the analysis. 



Cosmological parameters from RCS 



20<Rc<24 20<Rc<22 22<Rc<24 



FTTTl 1 I I I Mill 1 I I IIIW-IIII 1 I I I I Mill 1 I I I IIMI Mill 1 I I I Mill 1 I I I Mil 




2x10 



-5 _ 



10-5 




Is 




rriiil I I I 1 1 ml I I I 1 1 inrTiiil i i i 1 1 mi i i i 1 1 iTTI 1 1 nl I I I 1 1 ml I I I 1 1 III 
ffTTTj 1 I I I llllj 1 I I IIIHjllllj r 

2x10-5 




JilU I 



iiii 



rTTTTTl — I I I mill 



I I 



rnrr — i i i iiiiii — i i i iiii 



ILIU I I I Mini I 



1 10 

^[arcmin] 



10 



Fig. 1. — The upper panels show the measured variance of the aperture mass (M^p) as a function of aperture size for different samples of 
source galaxies. This signal corresponds to the "E'-mode. The lower panels show the variance (M^) when the phase of the shear is increased 
by 7r/2, and corresponds to the "B"-mode. The error bars indicate the 1<t statistical uncertainty in the measurements, and have been derived 
from the field-to- field variation of the 13 patches (thus the error bars include cosmic variance) . Note that the points are slightly correlated. 
We detect a significant "B"-mode on scales 5 — 10 arcminutes. On scales larger than 10 arcminutes the "B"-mode vanishes. The sample 
of bright galaxies (20 < Rq < 22) should not be affected significantly by systematics because their sizes are large compared to the PSF. 
Therefore the significant "B"-mode at scales of a few arcminutes is likely to be caused by intrinsic alingmcnts. To minimize the effect of 
intrinsic alignments on our cosmological parameter estimation, we will use the sample of galaxies with 22 < Rc < 24 to this end. 



We use single exposures in our analysis, and conse- 
quently cosmic rays have not been removed. However, 
cosmic rays are readily eliminated from the photometric 

catalogs: small, but very significant objects are likely to 
be cosmic rays, or artifact from the CCD. The peak finder 
gives fair estimates of the object size, and we remove all 
objects smaller than the size of the PSF. 

The objects in this cleaned catalog are then analysed, 
which yield estimates for the size, apparent magnitude, 
and shape parameters (polarization and polarizabilities) . 
The objects in this catalog are inspected by eye, in order 
to remove spurious detections. These objects have to be 
removed because their shape measurements arc affected 
by cosmetic defects (such as dead columns, bleeding stars, 
halos, diffraction spikes) or because the objects are likely 
to be part of a resolved galaxy (e.g., HIT regions). 

To measure the small, lensing induced distortions in the 
images of the faint galaxies it is important to accurately 
correct the shapes for observational effects, such as PSF 
anisotropy, seeing and camera shear; PSF anisotropy can 
mimic a cosmic shear signal, and a correction for the seeing 
is required to relate the measured shapes to the real lens- 



ing signal. To do so, we follow the procedure outlined in 
Hoekstra et al. (1998). We select a sample of moderately 
bright stars from our observations, and use these to char- 
acterize the PSF anisotropy and seeing. We fit a second 
order polynomial to the shape parameters of the selected 
stars for each chip. These results are used to correct the 
shapes of the galaxies for PSF anisotropy and seeing. 

The effect of the PSF is not the only observational dis- 
tortion that has to be corrected. The optics of the camera 
stretches the images of galaxies (i.e., it introduces a shear) 
because of the non-linear remapping of the sky onto the 
CCD. We have used observations of astrometric fields to 
find the mapping between the sky and the CCD pixel co- 
ordinate system, and derived the corresponding camera 
shear, which is subsequently subtracted from the galaxy 
ellipticity (see Hoekstra et al. 1998). 

As discussed in Hoekstra et al. (2002a) we use galaxies 
with 20 < Rc < 24 in our analysis. This selection yields a 
sample of 1,773,543 galaxies for which we could determine 
useful shape parameters. We note that the sample is not 
complete at the faint limit. 

The upper panels in Figure 1 show the measured vari- 
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ance of the aperture mass (M^p) as a function of aper- 
ture size 9. For all three samples we find that the 
"B"-mode (lower panels) vanishes on scales larger than 
~ 10 arcminutes, suggesting that neither observational 
distortions or intrinsic alignments of sources have cor- 
rupted our measurements. This is quite different from the 
VIRMOS-DESCART results presented by Van Wacrbcke 
et al. (2002) who found an almost constant "B"-mode out 
to scales of 30 arcminutes. 

To study possible residual systematics introduced by the 
correction for the PSF and camera shear, we split the data 
into various subsets: good/bad seeing, small/large PSF 
anisotropy, center/edge of chip. In all three cases, we find 
excellent agreement between the subsets. 

On smaller scales we detect a significant "B" mode. In- 
terestingly, we detect also a significant signal for the bright 
galaxies (20 < Rc < 22). These galaxies have sizes that 
are large compared to the PSF, and therefore they are less 
affected by residual systematics, whereas intrinsic align- 
ments are expected to be particularly important for bright 
galaxies. The selection of relatively bright galaxies in the 
Rc band results in a relatively narrow range in redshift 
as one essentially observes intrinsically fainter galaxies at 
redshifts z = 0.3 — 0.5. Only for fainter magnitude limits 
one overcomes the 4000A break, and more distant galax- 
ies enter in the sample. Hence, the "B"-mode at a few 
arcminutes is likely to be caused by intrinsic alignments. 
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Fig. 2. — The measured variance of the aperture mass (M^^) as 
a funetion of aperture size 9 for source galaxies with 22 < Rc < 24. 
The error bars have been increased to account for the unknown 
correction for the "B"-mode observed in Figure 1. 

The full sample (20 < Rc < 24) includes the bright 
galaxies, and hence the contribution of intrinsic alignments 
to the signal is larger compared to the faint (22 < Rc < 
24) sample. The fainter galaxies yield a slightly higher 
lensing signal, because the average redshift of the sources 
is higher. In addition the redshift distribution of galaxies 
with 22 < Rc < 24 is also broader (see, e.g., Hoekstra 



2001), which reduces the effect of intrinsic alignments. 

Currently it is not clear how to correct the "E"-modc, 
given the observed "B"-mode. We follow Van Waerbeke 
et al (2002) and add the "B"-mode in quadrature to the 
uncertainty in the "E"-mode, as this provides a conserva- 
tive limit to how well we can determine the lensing signal. 
The best signal-to-noise is obtained for the faint sample, 
and we will use this sample below for the cosmological 
parameter estimation. 

Figure 2 shows the resulting measurement of (Affp) as 
a function of scale when the error bars are increased using 
the observed "B"-mode. 

4. CONSTRAINTS ON COSMOLOGICAL PARAMETERS 

4.1. Source redshift distribution 

We now need to relate the observed lensing signal, pre- 
sented in Figure 2, to the cosmological parameters. To do 
so, we compute model predictions for a range of parame- 
ters using Eq. (8). However, the evaluation of Eq. (8) re- 
quires knowledge of the redshift distribution of the sources. 

Compared to other, deeper, cosmic shear studies, the 
RCS data have the major advantage that the redshift dis- 
tribution of the galaxies is known fairly well. Cohen et 
al. (2000) have determined spectroscopically the redshift 
distribution of galaxies down to a limiting magnitude of 
Rc = 24. However, their result is likely to suffer from 
sample variance. Comparison of the Cohen et al. (2000) 
redshifts with photometric redshifts (e.g., Fernandez- Soto 
et al. 1999) showed that the latter provide a suitable way 
to derive the redshift distribution of the sources we con- 
sider here. 

To do so, we have to account for the fact that the uncer- 
tainty in the shape measurements depend on the apparent 
magnitudes (and thus on the redshifts) of the sources: the 
contribution of distant, small faint galaxies (with noisy 
shape measurements) to the measured lensing signal is 
smaller compared to brighter galaxies. The relative weight 
as a function of magnitude (see Figure 7 in Hoekstra et al. 
2002a), is used to derive the "effective" redshift distribu- 
tion. It turns out that this weighting scheme changes the 
redshift distribution only slightly. We refer to Hoekstra 
et al. (2000; 2002a) for a more detailed discussion of the 
photometric redshift distribution used here (also see Van 
Waerbeke et al. 2002). 

It is expected that otir knowledge of the redshift dis- 
tribution of galaxies brighter than Rc — 24 will improve 
significantly in the near future. The current distribution is 
based on a relatively small sample of galaxies, and hence 
the data allow for a (small) range in redshift distributions. 

To allow for this uncertainty in the redshift distribution 
we marginalize our model predictions over the redshift dis- 
tribution allowed by the measurements of Fernandez-Soto 
et al. (1999). To do so, we parametrize the redshift dis- 
tribution using 



Pb{z) = 



exp 



(15) 



We fit this distribution to the observed "effective" red- 
shift distribution, fixing a = 4.7 and (3 = 1.7, which yields 
a best fit value of Zs = 0.302. We compute models for 
Zs G [0.274,0.337], which correspond to the ±3cr range 
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Fig. 3. — Contraints on and ag from RCS data only. The likelihood contours have been derived by comparing the measurements to 
CDM models with n = 1. We have marginalized the models over T and the Zg (for which we assumed a flat prior). The contours indicate 
the 68.3%, 95.4%, and 99.7% confidence limits on two paramaters jointly. Panel (a) shows the results marginalized over Ha; (b) gives the 
constraints for a flat cosmology (as suggested by the recent CMB results); (c) the results for an open (Qa = 0) cosmology. 



indicated by of the fit to the photometric redshift dis- 
tribution. For this clioicc of parameters the range in mean 
source redshift is (z) = 0.54 — 0.66. 

4.2. Parameter estimation 

In this section wc briefly discuss how we obtain con- 
straints on some of the cosmological parameters. The rel- 
evant parameters for the analysis presented here are the 
mass density Qm, the normalization of the power spec- 
trum ag and the shape parameter T (note that F ~ flmh 
in CDM cosmologies). 

We also consider Oa, although the RCS results are not 
expected to provide strong constraints on because 
lensing is not sensitive to this parameter. However, the 
combination of the RCS measurements with the VIRMOS- 
DESCART resuhs (Van Waerbeke et al. 200ia; 2002) 
can provide stronger constraints, because the two surveys 
probe the matter power spectrum at different redshifts. 
Hence, such a combined analysis enables a measurement 
of the growth of structure, and as a result it can break 
some of the degeneracies limiting our analysis. Further- 
more, we can impose priors on some of the parameters 
based on external measurements, such as the CMB or red- 
shift surveys. 

We compute model predictions on a grid, which are 
subsequently compared to the data. We have limited 
the model parameters to a realistic range. We consider 

n e [0,1], nA e [o,i], o-g g [0,2], and f e [0.05,0.5]. As 

mentioned above, we use Zg G [0.274, 0.337] to account for 
the uncertainty in the source redshift distribution. 

The measurements of (M^p) at different scales are 
slightly correlated. In order to compare the measurements 
to the models we have to compute the covariance matrix 
C which is given by 

Cij = {{di - iJ,i)'^{dj - iij)), (16) 

where di is the measurement at scale 9i, and /Xj is the 
"true" value for this scale (for which we take ensemble 
averaged value from the data). We compute the covari- 
ance matrix from the variance in the measurements for 



the 13 different patches, and hence it includes the cosmic 
variance. For a given model, we can compute the (log- 
likelihood) 

X'' = {di-mi)C-\di-mif, (17) 

which can be used to derive the likelihood for each set of 
parameters. From the measured values it is straight- 
forward to calculate the confidence contours. 

4.3. Joint constraints on flm and as 

We first consider the constraints we can derive on the 
combination of Cl„i and erg from the RCS data alone. With- 
out external priors these two cosmological parameters are 
degenerate. However, studies of cluster abundances give 
similar constraints, and it is therefore useful to compare 
the lensing results to the latter, popular technique. 

Figure 3 shows the joint constraints on fi,,,, and erg- We 
have marginalized over F and Zs, for which we assumed 
flat priors (i.e., all values have equal likelihood), which is 
rather conservative. Panel a shows the most general result, 
where we marginalize over f^A as well. For this choice of 
priors, we obtain erg = OA5to!lli^~°-^^ (95% confidence). 

The current CMB measurements favor a flat cosmology 
(e.g., de Bernardis et al. 2000; Jaffe et al. 2001; Netterfleld 
et al. 2001; Pryke et al. 2002; Stompor et al. 2001), and 
the results for these models are presented in panel b. For 
reference we also show the results for an open = 0) 
universe. 

The current constraints on flm and cis in Figure 3 are 
degenerate from lensing alone. The degeneracies can be 
broken by measuring the cosmic shear signal in both the 
linear and the non-linear regime (Jain & Seljak, 1997). 
However, the accuracy of our measurements at small scales 
(non- linear regime) is limited by the observed "B"-mode 
in the data. Hence, due to the increased uncertainties on 
small scales, we cannot break the degeneracy. 

However, if the "B"-mode is caused by intrinsic align- 
ments of galaxies, forthcoming multi-color data can be 
used to derive crude photometric redshifts, and we might 
be able to remove the contribution of intrinsic alignments 
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by correlating the shapes of galaxies at different redshifts. 
In addition, the planned CFHT Legacy Survey, which is a 
170 deg^ deep, multi-color survey will allow us to measure 
the cosmic shear signal with unprecedented accuracy, and 
these data will allow a determination of F, erg, and Qa 
from the shear two-point functions alone. 

We also note that in order to break the degeneracies, and 
infer the correct values for the various cosmological param- 
eters, the prescription for the evolution of the non-linear 
power spectrum needs to be improved. Van Waerbeke et 
al. (2001b, 2002) showed that the non-linear predictions 
fail for some cosmologies. The current prescription (Pea- 
cock & Dodds 1996) is not accurate enough for the future 
measurements. 

Fortunately, cosmological parameters can be estimated 
through a wide range of techniques, such as the fluctua- 
tions in the CMB and galaxy rcdsliift surveys. We can 
use the most recent results as priors to obtain tighter con- 
straints. 

In Figure 4 we show the constraints on Qm and (Tg, 
when we use Gaussian priors for the source redshift distri- 
bution Zs, F as determined from the 2dF galaxy redshift 
survey (Peacock et al. 2001; Efstathiou et al. 2001), and 
ritot = + from the Boomerang CMB measurements 
(Netterfield ct al. 2001). The contours have tightened con- 
siderably, and the data favor a low value for fi^- For this 
choice of priors, we find trg = (0.46tg:gf )0-°-52 (95% con- 
fidence) . 

We have also investigated whether the RCS data can be 
used to constrain F. Using our set of Gaussian priors, we 
find that we can only place a lower bound on F. We find 
a 95% confidence lower limit of F > 0.1 -|- 0.16flm- 



0-8 



1.5 



0.5 




Fig. 4. — Contraints on Qm and as using Gaussian priors on T, 
^tot = f2m -|- fiAl and the source redshift distribution Zs (which 
gives (z) = 0.59 ± 0.02). The likelihood contours have been derived 
by comparing the measurements to CDM models with Us = 1. For 
r we used the constraints from the 2dF survey F = 0.21 it 0.03 (Pea- 
cock et al. 2001; Efstathiou ct al. 2001), and for r2tot wc used the 
Boomerang constraints Qtot = 1.02 it 0.06 (Netterfield ct al. 2001). 
The contours indicate the 68.3%, 95.4%, and 99.7% confidence lim- 
its on two paramaters jointly. 



0.94t°o^ 



priors 



flat prior on Zs 

flat prior on Zg, flm -|- f^A = 1 

flat prior on z^, Q\ = 
Gaussian priors on Zs, F, and fitot 



Table 1 

Values of as and 68% confidence intervals for difFerent priors, 
assuming Qm = 0.3. 



(a) 


Cosmic shear survey 






RCS 

Bacon et al. (2002) 
Refregier et al. (2002) 


o.97;°:j° 

0.94 ±0.14 




Van Waerbeke et al. (2002) 


0.98 ±0.06 


(b) 


Cluster abundance 


0-8 



Carlberg et al. (1997) 
Bahcall & Fan (1998) 
Pen (1998) 
Borgani et al. (2001) 
PierpaoU et al. (2001) 
Reiprich & Bohringer (2001) 
Seljak et al. (2001) 
Viana et al. (2001) 



0.93 ± 
0.82 ± 
1.18t° 
1.00 ± 
0.72i° 

imto 

0.75 ± 
0.61 ± 



0.07 
0.03 

24 
22 

0.09 

07 
.05 
.07 
08 
08 
.06 

0.06 
0.05 



Table 2 

(a) Values of eg and 68% confidence intervals as derived from 4 
independent cosmic shear measurements (adopting i7m = 0.3, 
Qa = O-^, and F = 0.21; (b) determinations of as using cluster 
abundances. 



For comparison to other work, we list the 68% confi- 
dence limits on erg for = 0.3 for the different priors 
from our work in Table 1. It is also interesting to compare 
our results to other measurements. In Table 2a wc list our 
measurements with the most recent measurements of erg 
from weak Icrisiiig by large scale structure. The agreement 
between the four independent results is remarkable, in par- 
ticular given the fact that the measurements are based on 
a wide variety of data sets, using different telescopes, fil- 
ters, and depths. We can compute the ensemble averaged 
value of cTg from weak lensing, for which we find a value of 
0.92 ±0.03. 

We note that the value of CTg determined from the RCS 
data is essentially constrained by the mcasuremeiits at 
scales larger than 10 arcminutes, where the "B"-mode is 
negligible. The result from Van Waerbeke et al. (2002) 
might be biased high, because of their large scale "B"- 
mode. Both Bacon ct al. (2002) and Refregier et al. (2002) 
do not separate their signal into "E" and "B" -modes, and 
therefore it might also include some residual systematics. 

A widely used method to determine the normalization 
of the power spectrum uses the number density of rich 
clusters of galaxies (e.g., Borgani et al. 2001; Carlberg 
ct al. 1997, Eke, Cole & Frenk 1998; Bahcall & Fan 
1998; Pen 1998; PierpaoU, Scott, & White 2001; Reiprich 
& Bohringer 2001; Seljak 2001; Viana, Nichol, & Liddle 
2001). Such systems are rare, and as a result a very sen- 
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sitive probe of its, provided one can determine their mass. 
Tabic lb lists the results from 9 of these studies. The de- 
rived values for as from this technique range from values 
as low as 0.61 ± 0.05 (Viana et al. 2001) to values aroimd 
unity (e.g., Bahcall & Fan 1998; Pen 1998; Pierpoali 2001). 
Unfortunately, it is currently unclear why some of the re- 
cent values from cluster abundances are this low. From 
Table lb it is clear that there is a wide spread in values, 
whereas the statistical error bars are small. This points to 
underlying systematics, which complicate the determina- 
tion of as from cluster abundances. 




T 



reionization redshift. 

The CMB results arc taken from the compilation by 
Wang, Tegmark, & Zaldarriaga (2001), which is a combi- 
nation of many experiments, such as CBI (Padin et al. 
2001), BOOMERANG (Netterfield et al. 2001), DASI 
(Halverson et al. 2002), and MAXIMA (Lee et al. 2001). 
We use these results to obtain joint constraints on r, 
and as, adopting a baryon density u)b = 0.02, h = 0.7, 
f2„ = 0.3, and Qa = 0.7. We use the CMBFAST code 
(Seljak & Zaldarriaga 1996) to compute the theoretical 
CMB power spectra, and use the window function and 
covariance matrix provided by Wang et al. (2001). 

In Figure 5 we present the joint constraints on r and 
ris, where we marginalized over the value of as, using the 
prior as = 0.92 ± 0.03, as determined from weak lensing. 
The other cosmological parameters were fixed. The dot- 
ted contours correspond to the constraints derived from 
the RCS measurement of as alone. We stress that the 
results depend on the assumed priors (many cosmologi- 
cal parameters are fixed), and a full combined analysis of 
weak lensing and CMB data, marginalizing over all param- 
eters, is required for definite results. However, the results 
presented here demonstrate that such a combined analysis 
can be fruitful. 

For our choice of parameters, which arc close to those 
favored by the vast amount of data currently available, we 
can derive constraints on Ug and r. Marginalizing over 
r, we find = 1.02tQ Qg with 95% confidence. For r we 
obtain t = 0.12 ±0.09 (95% confidence, marginalized over 
ns). 

The reionization optical depth can be related to the red- 
shift of reionization. To do so, we use (e.g., Hu 1995; Grif- 
fiths & Liddle 2001) 



Fig. 5. — Joint constraints on the reionization optical depth r and 
the scalar spectral index Us from the combination of CMB results 
(compilation by Wang et al. 2001). We have marginalized over the 
value of (78, using a prior as = 0.92 it 0.03 as determined by the 
4 most recent weak lensing measurements. The other cosmological 
parameters were fixed as indicated. The dotted contours indicate 
the constraints from RCS alone (erg = 0.86 ± 0.05). The contours 
correspond to the 68.3%, 95.4% and 99.7% confidence limits on the 
two parameters jointly. 



4.4. Limits on the reionization optical depth t 

The measurements of the fluctuations in the CMB also 
give joint constraints on Q.^ and as- Recently, several 
papers argue that the CMB results yield a low normal- 
ization for the matter power spectrum (e.g., Lahav et al. 
2001; Melchiorri & Silk 2002), with values ranging from 
0.73 - 0.8. 

However, the current CMB data give relatively weak 
constraints on as, because the amplitude of the CMB 
power spectrum also depends on the reionization optical 
depth r, and the scalar spectral index Ug. Once the uni- 
verse is reionized, CMB photons can be Thomson scattered 
off of a free electron. As a result, the accousic peaks in the 
CMB power spectrum are suppressed. Although the effect 
is scale dependent (it does not affect the largest scales), 
it introduces a degeneracy between as and r. The matter 
power spectrum probed by weak lensing is not affected, 
and hence, the combination of the CMB and weak lensing 
results, in principle, can constrain r, and consequently the 



r(^) 



4.61 X 10" 



1 - 5 I X, 



2 / n 



(18) 



which is valid for flat cosmologies. Here Yp = 0.24 is the 
primordial mass fraction in Helium, and Xg is the ionized 
fraction. For a fully reionized universe {Xe = 1) the resul- 
tant r as a function of reionization redshift is presented in 
Figure 6. We have also indicated the 68% and 95% con- 
fidence limits on r (hatched regions), based on our con- 
straints on r. This allows us to read off the reionization 
redshift ^reion, for which we find Zreion — 14^9 with 95% 
confidence. 



5. CONCLUSIONS 

We have analysed ~ 53 square degrees of i?c-band imag- 
ing data from the Red-Sequence Cluster Survey (RCS), 
and measured the excess correlations in the shapes of 
galaxies on scales out to one degree. To this end we tise the 
aperture mass statistic, which allows a natural separation 
of the signal in a curl-free "E"-mode (which is expected 
from gravitational lensing) and a "B"-mode which pro- 
vides an important measure of residual systematics and 
the contribution from intrinsic alignments. 
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0.3 



Fig. 6. — The reionization redshift as a function of optical depth r 
as given by Eq. 18 (solid line) . The hatched regions indicate the 68% 
and 95% confidence limits on t. This corresponds to a reionization 
redshift of z,.; 



14^^ (95% confidence). 



For the lensing analysis we have selected a sample of 
galaxies with 20 < Rc < 24. On scales larger than 10 ar- 
cniinutes the "B"-niode vanishes, and therefore we are able 
to measure the cosmic shear signal with high accuracy on 
these scales. On smaller scales we find a small "B"-mode. 
We split the sample in bright (20 < Rc < 22) and faint 
(22 < Rc < 24) galaxies, and find a significant "B"-mode 
for the bright sample. These galaxies are large compared 
to the PSF, and therefore are not very sensitive to errors 
in the corrections for observational distortions. However, 
these galaxies probe only a relatively small range in red- 
shift, and intrinsic alignments of galaxies are expected to 
be important. Intrinsic alignments can introduce a "B"- 
mode and we conclude that the observed signal is likely to 
be caused by this efi'ect. A multi-color catalog, which is 
currently being created, will provide stronger constraints 
as it provides a better redshift discrimination of the galax- 
ies. 

The faint sample also shows a small "B"-mode, but it 

is smaller than the "E"-mode signal. It is not clear what 
is the cause of this signal, and hence the correction is un- 
certain. As a conservative approach we add the amplitude 
of the "B"-mode signal to the uncertainty in the measure- 
ment of the lensing signal. These measurements are then 
used to derive constraints on cosmological parameters. 
The relevant cosmological parameters for the analysis 



presented here are the mass density Qm, the cosmological 
constant Q\, the normalization of the power spectrum as 
and the shape parameter F (note that F ~ flmh in CDM 
cosmologies). To relate the model predictions to the data 
we use photometric redshift distributions for our sample 
of source galaxies (22 < Rc < 24), and marginalize over 
the uncertainty in the source redshifts. 

We investigated whether the RCS data can be used to 
constrain F. Using Gaussian priors for F (from 2dF) and 
Q,m + (from CMB), and the redshift distribution, we 
find that we can only place a lower bound on F for which 
we find F > 0.1 + 0.16f^,„ (95% confidence). 

We derive joint constraints on and as, by marginal- 
izing over F, Qa and the source redshift distribution, using 
different priors. Marginalizing over F and Q^, and using 
a flat prior for the source redshift distribution, yields a 
conservative constraint of as = 0A5t.Q'^ilil^'^° (95% con- 
fidence). A better constraint is derived when we use the 
Gaussian priors, which yields as ^ 0.46t|5:o?f^™°-^^ (95% 
confidence) . 

Comparison of the RCS results with three other re- 
cent cosmic shear measurements shows excellent agree- 
ment (Bacon et al. 2002; Rcfregier et al. 2002; van Waer- 
beke et al. 2002), and we find an ensemble averaged value 
of as = 0.92 ± 0.03 (68% confidence). The value of as 
determined from the RCS data appears to be free of sys- 
tematics. It is not clear, however, whether the results from 
the other surveys are biased high as a result of residual sys- 
tematics. Some recent studies (Reiprich et al. 2001; Seljak 
et al. 2001; Viana et al. 2001) of the abundance of rich 
clusters suggest low values for as, which disagree with the 
lensing results with high confidence. The reason for this 
discrepancy is currently unclear. 

The weak lensing results are also in good agreement 
with CMB measurements, when we allow the reionization 
optical depth r and the spectral index to vary. We 
have presented a simple demonstration of how the weak 
lensing results can be used as a prior for the parame- 
ter estimation from CMB measurement. We found that 
we can derive constraints on r and Ug (fixing flrn = 0.3, 
J7a = 0.7, ujb = 0.02, and h = 0.7). Doing so, we obtain 
Us = l-02to% and r = 0.12 ± 0.09, both with 95% con- 
fidence. We can relate the value of r to the reionization 
redshift ^reion for which we find .Zreion = 14lg with 95% 
confidence. We stress, however, that the results depend 
on the assumed priors, and that a full combined analy- 
sis of weak lensing and CMB data, marginalizing over all 
parameters, is required for definite results. Our prelimi- 
nary results, however, are not inconsistent with the latest 
observations of quasars at 2: ~ 6, which indicate that the 
universe may be just exiting the reionization epoch at that 
time (Barkana 2002; Becker et al. 2001; Djorgovski et al. 
2001; Fan et al. 2002). 
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